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Attractor solutions that give dynamical reasons for dark energy to act like the cosmological con- 
stant, or behavior close to it, are interesting possibilities to explain cosmic acceleration. Coupling 
the scalar field to matter or to gravity enlarges the dynamical behavior; we consider both couplings 
together, which can ameliorate some problems for each individually. Such theories have also been 
proposed in a Higgs-like fashion to induce gravity and unify dark energy and dark matter origins. 
We explore restrictions on such theories due to their dynamical behavior compared to observa- 
tions of the cosmic expansion. Quartic potentials in particular have viable stability properties and 
asymptotically approach general relativity. 



I. INTRODUCTION 

Over a decade ago two supernovae groups, the Su- 
pernova Cosmology Project and the High-Z Supernovae 
Search Team, provided evidence for an accelerated ex- 
pansion of the Universe [l[ . In recent years this discovery 
has gained more evidence from a variety of observations: 
further supernova data 0, measurements of the cosmic 
microwave background radiation Q and galaxy surveys 
Q ■ One possibility to explain this acceleration is to intro- 
duce a new component within the dynamics of General 
Relativity (GR), either as a uniform constant (cosmologi- 
cal constant) or as a scalar field evolving along a potential 
(as in the inflationary scenario) . Since the current accel- 
eration seems to be a unique phenomenon at least 
since the time of primordial nucleosynthesis, this poses 
a fine tuning problem unless some dynamical attractor 
exists. 

A different possibility is to look for acceleration as aris- 
ing from modifications of gravity. In many cases this can 
be viewed as coupling scalar fields nonminimally to grav- 
ity, within the framework of scalar-tensor theories Q, 
an approach called extended quintessence 0]. Attractor 
mechanisms can work here as well 04111 , and also in 
the case where a scalar field is coupled to (dark) matter 

num. 

In the present work we investigate the influence of both 
couplings on the attractor dynamics. From a phenomeno- 
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logical point of view this enriches the phase space and 
also can help with problems that arise from one coupling 
or the other. By comparison with observations of the 
cosmic expansion behavior we can constrain the allowed 
parameter space. From a theoretical point of view sev- 
eral models can lend motivation to such a combination 
of couplings. 

Induced gravity [l7j is similar to standard scalar-tensor 
theories, but gravity is induced by a Higgs-like field. One 
motivation stems from Einstein's original ideas to incor- 
porate Mach's principle into GR, by which the mass of a 
particle should originate from the interaction with all the 
particles of the universe, and so the interaction should be 
the gravitational one since it couples to all particles, i.e. 
to their masses or energies. To realize a stronger rela- 
tionship with the material contents, Brans and Dicke @ 
introduced their scalar-tensor theory of gravity, making 
the gravitational coupling, that is Newton's constant, a 
scalar function determined by the distribution of the cos- 
mic content. 

On the other hand, in modern particle physics the 
inertial mass is generated by the interaction with the 
Higgs field; the successful Higgs mechanism also lies pre- 
cisely in the direction of Einstein's idea of produ cing 
mass by a gravitational interaction. One can show [18j 
that the Higgs field as source of the inertial mass of the 
elementary particles mediates a scalar gravitational in- 
teraction of Yukawa type between those particles that 
become massive as a consequence of the spontaneous 
symmetry breaking. Due to the equivalence principle, it 
seems natural to identify both approaches. For this rea- 
son, [lj| [2(| proposed a scalar-tensor theory of gravity 
where the Higgs field of elementary particles also plays 
simultaneously the role of a variable gravitational con- 
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stant, instead of the scalar field introduced by Brans and 
Dicke. 

Put another way, if there is an interaction between a 
scalar (Higgs) field and matter, and of course there exists 
coupling between matter and gravity, why not close the 
loop by incorporating a (non-minimal) coupling between 
the scalar field and gravity? Or conversely, if one explores 
scalar-tensor theories, and gravity couples with matter, 
why not include an explicit interaction between the scalar 
field and matter? 

In this work we employ these ideas merely to motivate 
the couplings. For example, one might try to identify the 
dark energy (DE) with a Higgs-type (but not the Higgs) 
field that is coupled to some dark fermion sector. Ac- 
cordingly, in our model, gravity, i.e. the Ricci scalar R, 
couples to a scalar field </> through the non-minimal cou- 
pling (f> 2 R. As the field evolves to its energy minimum the 
Higgs coupling might give rise to the mass of some dark 
fermion that would account for the dark matter (DM) of 
the model. The scalar field evolves to a constant, to gen- 
erate the mass, and simultaneously generates Newton's 
constant through the non-minimal coupling to gravity. 
The resulting theory is an induced gravity in which GR 
is dynamically obtained through a Higgs mechanism from 
a scalar-tensor theory [2(J ■ 

The proposed Higgs mechanism is at present hypothet- 
ical, but phenomenologically interesting since it can ac- 
count for the mass of the DM and since its field can 
act as a DE to accelerate the cosmological expansion. 
While many DM-DE interaction models aiming to unify 
the two quantities have difficulties in getting them si- 
multaneously to exist and match observations, this can 
sometimes be made easier by adding a third element, such 
as inflation [Il|,[22| or, as in the present work, gravity. In 
any case, apart from the motivation, exploration of the 
dynamics of the matter- and gravity coupled system is of 
interest. 

We begin by describing the general field equations 
in Sec. [Hi identifying the contributions of the gravity- 
scalar and scalar-matter couplings. This is then exam- 
ined for the homogeneous and isotropic FRW universe 
background in Sec. IIII1 including the evolution equations 
for the scalar field and matter. Section ITVl discusses the 
effective potential, illustrating it for a symmetry break- 
ing form. In Sec. [V] we solve for the cosmic evolution of 
the field and matter, for general classes of potential and 
coupling, leading to constraints on the allowed parameter 
space from the cosmic expansion behavior. 



II. FIELD EQUATIONS 

The scalar-tensor theory Lagrangian used here is simi- 
lar to the one studied in the past for inflationary dynam- 
ics M: 
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where <f> is a real scalar field and a is a dimensionless pa- 
rameter. We use the metric signature (H ). Rewrit- 
ing the Lagrangian as 
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one sees that, formally, the Einstein-Hilbcrt action with 
the standard Newton's constant Gjv corresponds to cb 2 —> 



GR 



l/(aGjv)- Note however that even if on average 
<j> = i'GR; the theory will still be distinct from GR because 
of perturbations in the field. 

Varying Eq. ((TJ), one obtains the field equations 
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For the dark matter (DM) La- 



where a semicolon stands for a covariant derivative and 
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grangian we begin with a general "Yukawa" coupling, 
£m = i> (ij^dfj, — /(</>)) tp, similar to in [H|,[l5[, and later 
explore some particular forms. The dark matter thus has 
a <fi dependent mass /(<jd). 

The scalar field fulfills the generalized Klein-Gordon 
equation 



deb 8tt <P ~ deb 



(4) 



Note the term involving aR arises from the non-minimal 
coupling to gravity, and the right hand side comes from 
the coupling to dark matter. 

Taking the trace of Eq. <j3j> and substituting it into 
Eq. (dJ) to remove R, one obtains 
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(5) 

The first term in the parentheses is the normal GR term 
for an uncoupled scalar field (although with an altered 
prefactor), the T+4V of the second and third terms stem 
from the nonminimal coupling to gravity (specifically the 
trace of the first line of Eq. ([3]) , and the last term is due 
to the DM-DE interaction. 

This equation can be recast as: 



(0) 



where a prime stands for partial derivative with respect 
to <b. The effective potential is 
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from which the effective mass of the Higgs particle will 
be identified later on. We write the argument tp explic- 
itly for Cm to remind that only the partial derivative 
with respect to the field is relevant. If V(<f>) = Vb</> 4 , i.e. 
there is no intrinsic mass, then one sees that the poten- 
tial terms cancel out and the effective mass is completely 
determined by the trace T from the nonminimal gravity 
coupling and the DM-DE interaction term Cm- 

It is convenient to define the energy momentum tensors 
associated with the different contributions to the right 
hand side of Eq. ©: 

= + ^<j>., u - \^ X ^ V (8) 
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where the first is the standard scalar field contribution 
and the second stems from the nonminimal coupling to 
gravity. Again, note that the GR limit in which 
vanishes is not a = but <f> — > const. Due to the cou- 
plings, the components of the energy momentum tensor 
are not individually conserved but the total energy mo- 
mentum tensor is: 
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= 0. 
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The conservation equation for matter can be derived from 
this (using the equation of motion (j4} for <fr and the Ein- 
stein equations (O), giving 
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do 



(11) 



We will use this equation in the following to determine 
the evolution of the matter component. 



III. FRW COSMOLOGY EQUATIONS 

To study the background evolution of the universe we 
use the FRW metric 



ds 2 = dt 2 - a 2 {t) 



dr 2 
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and assume that the DM fluid is given by T^ v = p u^u v 
with the energy density p = nf(<j)), where n is the num- 
ber density, and = dx^/ds = 5^ the comoving 4- 
velocity. The energy momentum tensor behaves as a 
pressureless dust perfect fluid because we assume the 
dark matter to consist of non-relativistic fermions. The 
equivalence between the description in terms of the spinor 
field lagrangian as specified in the previous section on 
the one hand, and the "dust" description on the other 
hand is demonstrated in [23| for a standard linear cou- 
pling (/ oc </>). There, it is shown how a fermion field 
interacting with a scalar field <p, with action 



S = 



-giii/r/^dnip - y{4> - </>,)VV0, 
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where y and 4>* are constant, is equivalent to a model of 
a classical gas of pointlikc particles, with action 



S = ~H J y(0- 0*)dsi, 



(14) 



in the limiting situation where the fermions' de Broglic 
wavelengths are much smaller than the characteristic 
length scale of variation of the <j> field. Following the 
steps of that demonstration one can see that it is valid to 
replace the fact or y (4> — 0*) by an arbitrary function of 
</>, see also Ref. [24] . In our case we start with a fermion 
field with a coupling f((f>)^)tp. Then, a valid effective ac- 
tion is 

S=-J2j ftt)dsi = - j d A x V=9 f(4>)n(x), (15) 



where 
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This is exactly the action of a dust fluid with energy 
density p = f(<p)n. 

Thus, the matter Lagrangian is proportional to f(4>), 
i.e. Cm = —T = —p cx /(</>). From the gravity field 
Eq. ([3]) one obtains the generalized Friedmann equations 
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where H = a /a is the Hubble parameter. 

Rearranging the terms, the acceleration equation is 
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The redundant Klein-Gordon equation © becomes 
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where 

Kfr(0) 



3a 



V{<t>) + P {4>) 



dtp 



(21) 

Since the field is coupled to matter (either employing a 
Higgs or other mechanism), the source term V e g' involves 
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the density of the interacting DM fluid as well as the 
potential of the DE field. 

The conservation equation, Eq. (fTT|) . yields 
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so that the matter behaves the same way as in GR with 
a DM-DE interaction [l3| • This equation can be directly 
integrated to give 



(23) 



where no is the DM number density at present. This 
simply tells us that DM particle number is conserved and 
the change in energy per comoving volume is purely due 
to the varying mass /(</>). 



IV. EFFECTIVE POTENTIAL 

The nonminimal coupling to gravity, and the coupling 
to matter, in addition to adding symmetry to the rela- 
tions between the scalar field, matter, and gravity, also 
can create a nonzero vacuum expectation value (vev) - 
an effective cosmological constant - that can adiabati- 
cally evolve. The nonminimal gravity coupling can also 
reduce the driving term V e s (this coupling gives rise to 
the negative term in Eq. l2f p . slowing the field down. This 
slow roll can often alleviate instabilities in coupled mat- 
ter perturbations 
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We can examine these influences in terms of the effec- 
tive potential of the theory, which alters the bare scalar 
field potential through the coupling to dark matter and 
to gravity. Using the Yukawa coupling / cx </> in Eq. (f2"Tj) , 
the effective potential becomes 
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Note that the density terms in Eq. (|21j) - the interact- 
ing DM term C M and the trace term T stemming from 
the nonminimal coupling (j> 2 R - cancel out since p is lin- 
early proportional to (j>. If furthermore V(<fi) = Vo</> 4 , 
then the effective potential vanishes and the dark en- 
ergy field, even coupled, is masslcss and acquires infinite 
range. This is because in this case the theory has no 
explicit mass scale in it and is thus scale invariant. 

For the Yukawa coupling, the effective potential does 
not run with the density, so it differs from - and is ac- 
tually simpler than - what happens in the GR case with 
matter coupling. The effective mass is only determined 
by the potential terms. The field acts like quintessence 
with a V e ft(4>) given by the full Eq. flU}. 

For a Landau- Ginzburg symmetry breaking form for 
the Higgs potential, 



(25) 



where A a is a dimensionless constant, pV2 is the mass of 
the field at the potential minimum, and the Higgs ground 
state v, such that V(v) = 0, is given by 



6^ 
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In induced gravity, the Higgs potential V(4>) generates 
a time varying gravitational coupling (as in scalar-tensor 
theory) 



G(4>) 



(27) 



as <j) rolls from an initial state to its ground state and 
thus for a given field value determined by the potential, 
a needs to be chosen such that G(</>) ~ Gn- For example, 
we might choose a such that G(0 m i n ) = Gat. 

Relating this to the particle physics of the Higgs mech- 
anism, one has 
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where M Pl = 1/a/Gat « 1.2 x 10 19 GeV is the Planck 
mass, Mb the boson mass, and g a coupling constant. If 
we were to consider the standard model Higgs, one has 
M b = M w = 80 GeV for the W-boson and g = 0.18, 
therefore a « I0 33 . Such a value is huge and would 
not pass cosmological constraints on a discussed below. 
Therefore, we consider another Higgs-like particle with a 
much larger mass than that of the Higgs of the standard 
model of particle physics. The parameter a will need to 
be determined through cosmological observations. 

The effective potential for the Landau-Ginzburg plus 
Yukawa coupling case is 
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and has a minimum at <j) — v at all times, at which 
the mass of Higgs particle is m 2 H = V e s"{4> = v ) = 
A aV 2 /(3[l + 3a/(47r)]). 

We show the difference between the bare potential and 
the effective potentials in induced gravity and in GR with 
Yukawa coupling in Fig. [TJ Since the matter density con- 
tribution cancels out in induced gravity, the effective po- 
tential is time independent, while the GR case is not. The 
induced gravity case gives a broad, nearly flat minimum. 
When the field is very slowly rolling, the modifications 
to the right hand side of Eq. ([3]) are small. 

In this Ansatz the Higgs field, rolling along the effective 
potential, plays the role of DE and the particle masses 
produced by the Higgs mechanism provide the DM of the 
model. However, as we will see the dynamics does not 
generally favor such a form for the potential. 
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FIG. 1: The Landau-Ginzburg potential in the uncoupled, 
GR case (solid curve), in the matter coupled but minimally 
gravity coupled (GR) case for two values of the redshift z 
(dotted lines), and in the nonminimally gravity coupled in- 
duced gravity (IG) case (dashed curve). The matter coupling 
is taken to be of the Yukawa form. 



V. DYNAMICAL EQUATIONS 

We now examine the dynamics of the dark energy, the 
effective equation of state of the dark matter, and the 
overall expansion behavior. 



A. Variables and Equations of Motion 

The dynamics of the scalar field, non-minimally cou- 
pled to gravity and to dark matter, arc described by an 
autonomous system of equations which can be solved in 
a straightforward manner. One can define the following 
set of variables r2a| : 
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/ 6/f3/4> is a function of the field, with j3 



where k 
3a/(47r). 

Assuming spatial flatness, i.e. k = 0, the autonomous 
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where F = din f /din (j). The quantity C encodes the 
information on the matter coupling (in the first term in- 
volving /) and the non-minimal gravitational coupling. 

The fractional matter density comes directly from 
Eq. (JT7J) as 
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To evolve the field </> appearing in A or in d\nf/d\n(f>, 
one uses the auxiliary equation 
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We define the effective equation of state for a compo- 
nent in terms of n 2 pi (since the sum of these quantities 
is conserved, see Eg. [TO]) , by 
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This takes into account effective pressure terms due to 
the interactions. In particular, the matter equation of 
state gains an effective nonzero term, due to the non- 
minimal coupling to both the scalar field directly and to 
gravity through the <j) 2 R term. One has 

VP 

Wm = — 

One has w m = only when the field is frozen (x = 0), 
which restores GR as discussed in SeaHJ or when / ~ <f) 2 , 
which counteracts the k 2 term in Eq. (|3"T|) . 

The equation of state for the scalar field is more in- 
volved: 
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which depends on the nonminimal gravity coupling 
through (3, on the interacting DM-DE term (F), and on 
the potential (U = d In V/d In . 
The total equation of state 
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takes the same form with or without coupling to matter. 
The deceleration parameter q = (1 + 3w to t)/2 as usual. 




B. Critical points 

Investigating the critical (or fixed) points of the dy- 
namical system, we find the relation (for y c ^ 0) 
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We restrict to the cases of U and F being constant in 
this analysis; otherwise the equations arc transcendental. 
The dark energy density and equation of state are 



tt<p,c = l-^-(4-U)x c ^ 
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The fixed point x c is given by a quadratic equation. 
The first solution is 
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This is a stable fixed point over a wide range of parame- 
ters (see next section) and it is independent of the mat- 
ter coupling f(4>). Note that we obtain a future dc Sitter 
state (w = — 1) for V ~ cf> 2 or V ~ <fi 4 asymptotic behav- 
ior. The U = 4 case represents an attractor to GR, since 
in this case x c = = (f>. An x c — solution was also 
found in GR with no matter couplings H, HH HH , and 
with an exponential coupling [161. and also in induced 
gravity with no matter couplings jlll [27j . For 2 < U < 4 
the fixed point gives a phantom [w c < —1) attractor; a 
phantom solution was obtained in Ref. [28| for F = and 
a particular choice of U. The general attractor solution 
is illustrated in Fig. [2j Note that there is a wide range 
of U and a (and all of F, which does not enter) where 
the attractor gives acceleration, and in fact stays close to 
w = — 1. In all cases a — > makes u^.d — > — 1. 

The second solution (mostly a saddle point) does de- 



FIG. 2: The dark energy equation of state for the first critical 
point is plotted as a function of the power law index U of the 
potential and the value a of the gravity coupling (the matter 
coupling does not affect this critical point). The value w = — 1 
is indicated by the green mesh plane. Note the attractor gives 
acceleration, and w ~ —1, over a wide range of values. 



pend on /(</>) and is given by 
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(49) 
(50) 



(except when F 



2 

also, since for F = U this solution does not exist). A 
solution within GR for an ex pon ential potential and an 
exponential coupling function |l2[ also has a critical point 
depending on the equivalent of F and U. 

Note that x 2 , y 2 , and f2 rn can all be greater than unity. 
This is possible because of the negative term —Hcjxfi in 
the modified Friedmann equation (|17[) . which permits the 
effective dark energy density to go negative. This causes 
no physical problems within our Ansatz, as we demon- 
strate in Sec. IV Dl 

Three further mathematical solutions exist for y c 
0. These have x C 4 ;C 5 = \fp ± \J\ + (3, which are mostly 
unstable and inaccessible, but of more interest is 



x C 3 



y^(l--F) 
3 + /3(2 + F) 

(3(F-l)[5 + F + 2f3{2 + F)} 

[3 + /3(2 + F)] 2 
/?(F-l)(F-2) 
3[3 + /3(2 + F)] ' 
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This solution is important because y ss corresponds to 
k 2 V/ (3H 2 ) w 0, which holds at early times. We indeed 
find a metastable attractor to this behavior at high red- 
shift (basically, dx/dN w but y keeps growing). Note 
that this high rcdshift attractor is generally iiu = w m , 
i.e. a scaling solution. However, V (and y) arc not ac- 
tually zero, just small compared to the other densities, 
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so the case F = 1, which would give x c — 0, breaks the 
condition y c <C x c and instead here the dynamics forces 
W0 = — 1 and w m = 0. 

C. Stability 

Carrying out a linear stability analysis, we find that the 
stability of the critical points generally depends on the 
values of a (i.e. /3), U, and F, making analytic statements 
difficult. However, a reasonable rule of thumb is that for 
values of these variables not too positive or too negative, 
the first critical point is a stable node (both eigenvalues of 
the perturbation matrix negative) and the second critical 
point tends to be a saddle point (one eigenvalue positive, 
one negative). See [25[ for a general discussion of stability 
analysis and classification. 

We show the 3-dimensional eigenvalue surfaces for the 
first critical point in Fig. |3l for three different values for 
a. The region of stability exists for a broad interval of 
U and F, including "natural" values. As a decreases, 
the area of stability grows, as seen in the sequence of the 
three plots. For a < 0.1 the system is stable for any 
U e [-20, 20], and F e [-10, 10]. 

By contrast, the second critical point is a saddle point, 
as seen in Fig. [4] One eigenvalue is always positive and 
the other negative, independent of the chosen U and F. 
As a decreases, the eigenvalue sheets tend to pull further 
away from the zero plane. Since the qualitative nature of 
one positive eigenvalue, one negative eigenvalue does not 
change, we show two values of a together in the figure. 
Also note that the eigenvalues diverge for F = U, as 
expected from Eq. (|48[) . 

The third critical point is also a saddle point for 
U £ [-10,20] and F € [-10,10], but it has a stabil- 
ity window for a = 1 in the region U G [—10, —20] and 
F E [—1, —5]. As a — > the sheets flatten, one below and 
above the fi — plane, making this a saddle point for any 
U and F. Finally, the fourth and fifth fixed points are 
mostly unstable, with small regions of stability (for big 
|F| and \U\) for a > 0.5, but these solutions are anyway 
uninteresting since they are nonaccelerating and require 
y = = V while fi^ = 1, hence </> = and vanishing 
gravity. 

A key point is that U = 4 always has a stable attractor 
to the first critical point, with w = — 1 and = 0, for all 
values of a and F (even F ^ constant, as we discuss in 
the next section). So for a quartic potential there is an 
attractor leading to ACDM and a restoration of GR. 

D. Dynamical History 

Turning now to the full dynamics, some numerical so- 
lutions are illustrated in Fig. [5] At high redshift, where 
H 2 is dominated by the matter density, y is small and the 
field quickly forgets the initial conditions and goes to the 
scaling solution of the third critical point (except when 



a= 1.0 




FIG. 3: Surfaces for the two eigenvalues fii (upper red) and 
fi2 (lower blue) of the first critical point. The stability region 
for the critical point is where both eigenvalues are negative, 
i.e. below the (green horizontal) fi — plane, and is a function 
of the power law index U of the potential, F of the matter 
coupling, and the value a of the gravity coupling. The three 
panels shows how the surfaces change with a: the smaller a, 
the more stable the system. 
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a= 1.0,0.1 




FIG. 4: Surfaces for the two eigenvalues fii (upper red) and /j.2 
(lower blue) of the second critical point. Since one eigenvalue 
is positive, this critical point is not stable. As a decreases, 
the sheets pull away from the zero plane. 



F = 1, then = — 1). As the contribution of the dark 
energy potential energy becomes relatively more impor- 
tant, it then approaches the asymptotic attractor of the 
first critical point. 

We define the present (a = 1) by when fi^ = 0.72. For 
clarity we show in Fig. [5] the results for a value of the 
gravity coupling a = 1; for small a the deviation from 
w(a) = — 1 will scale roughly as a. The high redshift 
metastable attractor and the asymptotic future stable 
attractor behaviors can clearly be seen. Note that for 
the Yukawa coupling (F = 1) the dynamics stays close 
to w = —1. For F < 1, the —Hcjxj) term can drive the 
effective dark energy density through zero, making w go 
to ±00. This has no physical pathology, since the dark 
energy density is merely an effective quantity and as we 
will soon see the matter and total equations of state are 
well behaved. In the future, the attractor solutions for 
V ~ 4> 2 and V ~ t/> 4 , i.e. U = 2 and 4, are the de Sitter 
state w = —1. The latter case represents the attractor 
to GR. 

Examining the future attractor solution more closely, 
we see from Eq. ({47]) that within the range 2 < U < 4, 
the attractor is to a phantom state w < — 1, while out- 
side this range it is to w > —I (unless U gets extremely 
negative). In general, reasonable inverse power law po- 
tentials do not provide w w — 1 and would be disfavored 
by observations. Figure [6] shows the full behaviors for 
a variety of U values, fixing the matter coupling to the 
Yukawa form / ~ <fi (F = 1). 

Regarding the early time behavior, as stated the dark 
energy equation of state becomes undefined when the 
dark energy density passes through zero, as can occur for 
F < 1. Since the dark energy density vanishes, however, 
this has no physical effects. In particular, the matter 
equation of state w m and total equation of state u>tot ex- 
hibit no sign of this occurrence. In Fig. [7] we see that at 
high redshift w m goes to its metastable attractor solution 
(third critical point), and because of scaling w to t = w m . 
In the Yukawa coupling case (F = 1), w m = 0, as it is 
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FIG. 5: The dark energy equation of state is plotted for three 
power law indices of the potential and the matter coupling. 
All solutions asymptote to the stable fixed point w c i(U), in- 
dependent of F. Note both quadratic and quartic potentials 
(U = 2, 4) reach the de Sitter state w — —1. The deviations 
of w(a) from —1 scale linearly with a for small a. 
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FIG. 6: The evolution of the dark energy equation of state is 
shown for various power law potentials V ~ 4> u , with Yukawa 
coupling / ~ (f> and gravity coupling a = 1. The attractor 
values 1 + w c scale with a for small a. Note that 4> 2 an d </> 4 
potentials give w = — 1 attractors, while power law indices in 
between give phantom attractors. 
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for the F = 2 case always as well. Both w rn and Wtot 
then smoothly evolve toward the future stable attractor 
solution. Both w m and 1 + Wtot scale linearly with a for 
small a. For F = 1, all quantities are well behaved and 
follow ACDM until recently. 
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FIG. 7: The matter equation of state w m and total equation 
of state Wtot are plotted for the same potential and matter 
coupling cases as in Fig. [5] The asymptotic future values 
of w m and 1 + Wtot scale linearly with a for small a. For 
Yukawa coupling (F = 1), until recently w m was zero and 
wtot followed the standard ACDM history (solid black line, 
also for U — 4, F — 1) as well. 

When a is small, the equations of state of matter, dark 
energy, and the total energy are nearly the same as for 
ACDM, regardless of F and U. This is demonstrated in 
Fig. E for a = 0.1. 

A small value of a corresponds to <f> greater than the 
Planck energy to induce GV (see Eq. l2~T|) . However, there 
is no actual energy density that we are taking to be 
greater than the Planck scale, so it is not clear there is a 
problem in treating this as a low energy effective theory. 
In fact, our theory does not have an explicit Planck scale 
in it so it is not even clear what the cutoff energy is above 
which our theory is expected to break down. 

Another issue is how to apply standard scalar-tensor 
theory limits from the Solar system, e.g. on the Jordan- 
Brans-Dicke parameter lojbd- Without the matter cou- 
pling and with U = 0, one would say that lojbd = 2n/a 
and so these constraints [29[ require a < 1CP 4 . In Ref. 
[27j it is shown a compatibility of the model U = 4, 
F = with current cosmological parameters respecting 
Solar system constraints. However, the matter coupling 
might alter that conclusion. Ref. [30[ considers the differ- 
ences between matter coupled and nonminimally gravity 




a 



FIG. 8: As the U = 4 case of Fig. but for a = 0.1. Devi- 
ations of vim from zero, and u>tot from the ACDM behavior 
are less than 1CP 2 . 

coupled models, but so far no analysis has been carried 
out including both effects. The equations for the field 
perturbations with all the couplings are quite compli- 
cated; we leave that analysis for future work and here 
only consider effects on the cosmic expansion. 

We can consider going beyond constant F and U. In 
the case of an exponential matter coupling, the logarith- 
mic derivative F is not constant. Figure shows w(a) 
for such a case, where / ~ e so F = bcj> (here the field 
values are all in units of the Planck mass). The dynam- 
ics is now dependent on the initial value of the field <f>i, 
and resembles that of the constant F case with F = b<pi 
(a better approximation is evaluating this at a ^ part- 
way between c^i and a later time value). In particular 
this applies to the high rcdshift attractor. However, note 
that the late time attractor given by Eq. (TT7|) remains 
the same, independent of F and </>,. 

One can also consider potentials where the logarith- 
mic derivative U is not constant. Figure [10] illustrates 
w(a) for a potential of Landau-Ginzburg form (|25[) . In 
this case, with varying U, the dynamics is quite different 
from the constant U case, and no attractor solution is 
apparent. 



VI. CONCLUSIONS 

The dynamics of a scalar field model incorporating 
both nonminimal coupling to gravity and coupling to 
matter is rich. The relations to induced gravity through a 
Higgs mechanism and the symmetrical ideas of coupling 
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FIG. 9: Constant U potentials with exponential coupling / ~ 
e 6< *, shown for 6 = 2. For U = 4 (thick, black curves), the 
cases for three initial field values <f>i are shown; for (7 = 2 and 
U — —2 only the bcj>i = 2 case is plotted. The w(a) dynamics 
depends on <f>i (although it is quite similar to a corresponding 
constant F) but the final attractor depends on neither / nor 
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FIG. 10: Comparison of the Landau-Ginzburg (LG) form, 
with either constant F or non power law /, to the 4> 4 potential 
shows significant differences. In particular, the non power law 
potential does not exhibit attractor behavior. 



around the circle of scalar field, matter, and gravity are 
interesting. Furthermore, adding both couplings can re- 
lieve some problems of either individual coupling. 

A set of attractor solutions for the dynamics is found, 
with a stable attractor to a dark energy dominated uni- 
verse with w either equal to or near — 1. High redshift 
behavior is standard ACDM for a Yukawa matter cou- 
pling (better behaved than if there were no coupling), 
and a scaling solution otherwise. Throughout cosmic his- 
tory the matter and total equations of state can also be 
acceptably near ACDM behavior, for a not much larger 
than unity, depending on the couplings. 

A simple quartic potential V ~ <p 4 appears viable 
(again better behaved than if the potential were con- 
stant), with a stable de Sitter attractor for all values 
of gravity and matter couplings. Gravity is asymptoti- 
cally restored to general relativity (at least as far as the 
background behavior is concerned). A quadratic poten- 
tial also has a do Sitter attractor. 

We identified the roles of the various couplings in the 
evolution equations and the effective potential. The 
slowing of field dynamics due to the gravity coupling 
likely alleviates problems with coupled matter instabil- 
ities, through the adiabatic mechanism discussed by [lj|, 
and the coupling to matter may help issues with gravity 
tests, but the system of perturbation equations becomes 
quite complicated due to the two extra scalar field cou- 
plings and we leave that to future work. Here we concen- 
trated solely on the field and expansion dynamics, which 
in itself can constrain the parameter space, while point- 
ing the way to interesting attractor behaviors and the 
possibility of dynamically establishing w ~ — 1 . 



Acknowledgments 



We gratefully acknowledge a UC-MEXUS-CONACYT 
Visiting Fellowship for JLCC to spend a sabbatical at 
Berkeley JLCC thanks the Berkeley Center for Cosmo- 
logical Physics for hospitality and CONACYT for grant 
No. 84133-F. This work has been supported in part 
by the Director, Office of Science, Office of High En- 
ergy Physics, of the U.S. Department of Energy under 
Contract No. DE-AC02-05CH11231, and the World Class 
University grant R32-2009-000-10130-0 through the Na- 
tional Research Foundation, Ministry of Education, Sci- 
ence and Technology of Korea. 



11 



[1] A. G. Riess et al., Astron. J. 116, 1009 (1998); 

S. Perlmutter et al, ApJ, 517, 565 (1999). 
[2] R. Amanullah et al., ApJ, 716, 712 (2010). 
[3] E. Komatsu et al, arXiv: 1001.4538. 
[4] B. A. Reid et al., MNRAS 404, 60 (2010). 
[5] E.V. Linder and T.L. Smith, arXiv: 1009.3500. 
[6] C. Brans and R. Dicke, Phys. Rev. 124, 925 (1961). 
[7] F. Perrotta, C. Baccigalupi, and S. Matarrese, Phys. Rev. 

D 61, 023507 (2000). 
[8] N. Bartolo and M. Pietroni, Phys. Rev. D 61, 023518 

(2000). 

[9] C. Baccigalupi, S. Matarrese and F. Perrotta, Phys. Rev. 

D 62, 123510 (2000). 
[10] T. Chiba, Phys. Rev. D 64, 103503 (2001). 
[11] V. Faraoni and M.N. Jensen, Classical Quantum Gravity 

23, 3005 (2006). 
[12] L. Amendola, Phys. Rev. D 62, 043511 (2000). 
[13] S. Das, P.S. Corasaniti, and J. Khoury, Phys. Rev. D 73, 

083509 (2006). 

[14] C. Quercellini, M. Bruni, A. Balbi, and D. Pietrobon, 

Phys. Rev. D 78, 063527 (2008). 
[15] P.S. Corasaniti, Phys. Rev. D 78, 083538 (2008). 
[16] R. Bean, E.E. Flanagan, I. Laszlo, and M. Trodden, 

Phys. Rev. D 78, 123514 (2008). 
[17] P. Minkowski, Phys. Lett. B 71, 419 (1977) ; A. Zee, 

Phys. Rev. Lett. 42, 417 (1979); S.L. Adler, Rev. Mod. 

Phys. 54, 729 (1982) . 
[18] H. Dehnen, H. Frommert, and F. Ghaboussi, Int. J. Theo. 

Phys. 29, 537 (1990); H. Dehnen and H. Frommert, Int. 

J. Theo. Phys. 30, 985 (1991). 
[19] H. Dehnen, H. Frommert, and F. Ghaboussi, Int. J. Theo. 



Phys. 31, 109 (1992). 
[20] J. L. Cervantes-Cota and H. Dehnen, Phys. Rev. D 51, 
395 (1995); 

J. L. Cervantes-Cota and H. Dehnen, Nucl. Phys. B 442, 
391 (1995). 

[21] A. R. Liddle, L. A. Urena-Lopez, Phys. Rev. Lett 97, 
161301 (2006). 

[22] T. Padmanabhan and T. R. Choudhury, Phys. Rev. D66, 
081301(R) (2002); R. J. Scherrer, Phys. Rev. Lett. 93, 
011301 (2004); A. Arbey, Phys. Rev. D74, 043516 (2006); 
P. J. E. Peebles and A. Vilenkin, Phys. Rev. D59, 063505 
(1999); J. E. Lidsey, T. Matos, and L. A. Urena-Lopez, 
Phys. Rev. D66, 023514 (2002); T. Matos, J.-R. Luc- 
vano, and H. Garcia-Compean, hep-th/0511098. 

[23] G. R. Farrar and P. J. E. Peebles, ApJ, 604, 1 (2004). 

[24] T. Damour and A.M. Polyakov, Nucl. Phys. B 423, 532 
(1994); 

[25] E. J. Copeland, M. Sami, S. Tsujikawa, Int. J. Mod. 

Phys. D 15 1753 (2006). 
[26] T. Damour and K. Nordtvedt, Phys. Rev. D 48, 3436 

(1993); D.I. Santiago, D. Kalligas and R.V. Wagoner, 

Phys. Rev. D 58, 124005 (1998). 
[27] F. Finelli, A. Tronconi, and G. Venturi, Phys. Lett. B 

659, 466 (2008). 
[28] L. Perivolaropoulos, JCAP 10, 1 (2005). 
[29] B. Bertotti, L. less, and P. Tortora, Nature 425, 374 

(2003). 

[30] V. Pettorino and C. Baccigalupi, Phys. Rev. D 77, 
103003 (2008). 



